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Abstract 



o 

In this work we consider the evolution of the interactive dark fluids in the background 
of homogeneous and isotropic FRW model of the universe. The dark fluids consist of 
a warm dark matter and a dark energy and both are described as perfect fluid with 
barotropic equation of state. The dark species interact non-gravitationally through an 
additional term in the energy conservation equations. An autonomous system is formed 
Q-i in the energy density spaces and fixed points are analyzed. A general expression for the 

deceleration parameter has been obtained and it is possible to have more than one zero 
bJQ' of the deceleration parameter. Finally, vanishing of the deceleration parameter has been 

examined with some examples. 
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1 Introduction 

Recent observational data coming from Supernove la (SNe la) [1], large scale structure formation (LSS) 
[2], cosmic microwave background radiation (CMBR) [3], baryon acoustic oscillation (BAO) [4] and 
weak lensing [5] put a big question mark to the standard cosmology. As the gravity of both baryonic 
(ordinary) matter and radiation is attractive in nature, so standard cosmology predicts a deceleration 
of the universe while observationally the universe is currently going through an accelerated expansion. 
To overcome this challenging and intriguing problem in cosmology people are trying to modify the 
standard cosmology in two ways: (I) by incorporating a new form of matter called dark energy (DE) 
which has marginally dominant negative pressure component that drives the present acceleration of 
the universe [6]. (II) Secondly, the doubt about the applicibality of the General Relativity to describe 
the universe as a whole, results alternative theories of gravity. Also alternative theories of gravity are 
proposed to resolve the puzzles of theoretical and experimental gravity [7] . These include f(R) gravity 
theory [8], theories with other curvature invariants [9], Coupling the Ricci scalar to a scalar field [10] 
by introducing a vector field contribution [11], or by gravity theory in higher dimensional space-times 
[12]. In these theories to produce cosmic acceleration one does not require any exotic matter (i.e, DE) 
rather it is provided by the extra terms in the usual Friedmann equations on large scales [13]. 
Due to lack of any definite conclusion, various cosmological models have been proposed to match 
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with recent observed data. ACDM model is the simplest one in this series. Although the model 
predicts cosmic acceleration as well as a reasonable agreement with observational data, but there 
are some embarrassing issues related to this model namely cosmological constant problem [14] (the 
huge discrepency between the observed value of the cosmological constant and the one predicted in 
quantum field theory), coincidence problem [15] (although generically small, but the cosmological 
constant happens to be exactly of the value required to become dominant at the present epoch) and 
recently, it was shown that the ACDM model may also suffer the age problem [16]. 
To overcome these puzzles, it is natural to persue alternative possibilities to explain the mystery of 
dark energy. In the literature it is assumed that the energy of the vacuum is zero (by some unknown 
cancellation mechanism) and in its place a dark energy components having dynamically variable 
equation of state is considered. Over the last decade, various dark energy models have been proposed 
namely, Quintessence [17], Phantom [18], K-esssence [19], Tachyon [20], Quintom [21], chaplygin gas 
[22], Holographic DE [23], the new age graphic [24], the Ricci DE [25] and so on. Motivated by the 
supersymmetric field theories and String/ M theory these DE models are represented by some effective 
scalar fields. 

Further, though both the dark components (DE and DM) evolve differently, but cosmic observations 
predict their energy densities of the same order today. To overcome this Coincidence problem [26], 
both the dark components are assumed to interact non-gravitationally [27] through an additional term 
in the conservation equations. Also in the perspective of the present data the DE should be chosen 
such that there is a smooth transition across the phantom barrier from the above in near past [28]. 
In the present work, we propose a phenomenological cosmological model where the two dark species 
interact among themselves. The equation of state parameter for warm dark matter is chosen to be 
positive but gradually decreases while the DE equation of state parameter is negative throughout 
the evolution but otherwise unrestricted. The model describes the evolution of the universe from an 
initial dense radiation era (or a barotropic perfect fluid) when DM dominates over DE. Subsequently, 
DE takes an upper hand and the universe gradually enters into the accelerating phase. We mainly 
focus our attention how the evolution of the universe and its cosmic acceleration depend on the 
evolving equation of state parameters of the two dark specieses. we shall also examine (by examples) 
whether the accelerated expansion will continue forever or it is possible to have an ever expanding 
DE-dominated universe which passes through from deceleration to acceleration and vice versa. 
The paper is organized as follows: Basic equations for the phenomenological model has been presented 
in section 2. For two specific choices of the interaction term, the energy conservation equations are 
reduced to an autonomous system and critical points are analyzed in section 3. Section 4 presented 
a general prescription for the deceleration parameter and in section 5, it has been examined whether 
more than one transition is possible or not by some examples. Finally, in section 6, there is a brief 
discussion about the results obtained in this work. 



2 Basic equations: 

We start with homogeneous and isotropic FRW model of the universe filled with dark fluids- a warm 
dark matter of energy density p m and a dark energy component described by the dark energy density 
p d . As we are interested with the late time cosmic evolution so we have neglected other matter 
components namely radiation, baryons etc. For simplicity, both the dark fluids are assumed to have 
barotropic equation of state namely, 

Pm = UmPm and p d = UJ d p d (1) 

where u) m (> 0) and u d (< 0) are in general chosen as variables. So the Einstein field equations take 
the form 

3(H 2 + 4) = Pm + Pd (2) 
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and 

H-^ = --[ Pm {l + u d ) + Pd {l+u d )\ (3) 

(X Z, 

If we assume that the two dark components interact among themselves, then the continuity equations 
are of the form 

p m + 3Hp m (l + U} m ) = Q (4) 

and 

p d + 3H Pd (l + uj d ) = -Q (5) 

where Q stands for the interaction between the dark sectors. For the time being we do not use any 
specific choice for Q, only assume that Q does not change its sign during the cosmic evolution. Now 
the above continuity equations can be written in the non-interaction form with effective equation of 
state parameters as 

u me = u m - — — and uj de = u d + — — (6) 

3Hp m 3Hp d 

As usual, if we assume that the energy is transferred from dark energy (DE) to dark matter (DM), 
i.e, Q > 0, then we have the following observations: 

(a): For cold DM (i.e, uj m = 0) we have effectively some kind of exotic DM (i.e, -ve equation of state 
parameter) in the presence of interaction, (b): For warm DM (i.e, oj m > 0) it is likely to have a 
possible change of sign of the effective equation of state during cosmic evolution depending on the 
nature of the interaction and the strength of the coupling constant appeaaring in Q. (c) If the DE is 
chosen as the cosmological constant (i.e, co d = —1) then the effective DE will behave as a quintessence 
field (i.e, uj de > —1). Also the effective DE fluid may be quintessence era even we choose phantom DE 
(i.e, uj d < -1). 

We now introduce the dimensionless parameter r = ^ 2i , which is also termed as the coincidence 
parameter. Using the energy conservation equations (4) and (5), the evolution equation for r with 
respect tor = 31na is given by (after some simplification) 

TT =r[^-^) + - { -^L] (7) 

where pt = p m + p d , is the total energy density of the dark fluids. Also using the above energy 
conservation equations the evolution of pt is given by 



The fixed points of equations (7) and (8) are given by 

1 +u d 
rc = 

1 + W m 

and 



(9) 



Q LO m - LO d 

pc ~ m{i + uj m ){i + uj d ) ( ' 

where in the last equation we have used equation (9) for r c . It is to be noted that as r c and p c are 
positive definite so we must have phantom DE, i.e, uj d < — l.If oj m = 0, i.e, for cold DM to have a 
realistic fixed point (i.e, +ve stationary energy density and coincidence parameter) there must be an 
exchange of energy from phantom DE to cold DM. The same is true for warm DM also. It is to be 
noted that the result still remains valid even for — 1 < uj m < 0. 
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3 Fixed points and stability analysis for Holographic 
dark energy model: 

In this section we shall choose DE as the Holographic dark energy (HDE) model with event horizon as 
IR cut off. If Re is the radius of the event horizon then from the Holographic principle [29] we have 

where the dimensionless parameter V is determined by observational data [30]. However, in the 
present work we have taken c to be arbitrary. Also £l d = is the usual density parameter for DE. 
Now using the conservation equation (5) and the Friedman equations (2) and (3) the evolution of the 
DE can be written as 

— = (i - n d )[3uj m + h \--ft~ ( 12 ) 

U d c Hp d 

where / stands for differentiation with respect to x = lna.The equation of state parameter for the HDE 
has the form 

Note that oj d does not depend explicitly on uj m . Also it should be mentioned that if co m = then with 
proper choice of Q, one gets back to equations (5) and (6) of reference [30] (or equations (23) and (24) 
of reference [31]) from equations (12) and (13) respectively. We shall analyze the evolution equations 
for the following two choices of the interaction term : 

(a) :Q = 3Hb 2 {p m + Pd ) , (b) :Q = 3\Hp d (14) 



with b 2 and A as the couplng constants. 
Case I: Q = 3Hb 2 (p m + p d ) 

This expression for the interaction term was first introduced to study coupling between a quintessence 
scalar field and a pressureless cold DM field [27], so that the coincidence problem can be resolved by 
scaling solution and the universe approaches a stationary stage (i.e, r = constant). Subsequently, this 
form of interaction was derived in the context of HDE with Hubble scale as the IR cut off [27] and 
also with event horizon as the IR cut off [31]. Now for this choice of Q and using the equation of state 
parameter ui d for HDE from equation (13), the energy conservation equations can be written as 

P'm = Sp m [b 2 (l + -) - (1 + U m )] (15) 

and 

pi d = -2p d [l - —±=] (16) 
cyl + r 

If oj m is assumed to be constant then equations (15) and (16) form an autonomous system in the (p d , 
p m ) plane. Hence at the fixed points the coincidence parameter has the constant value 

r = \-l (17) 
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which implies c 2 < 1. Also uj m corresponding to the fixed points has the value 

Wm = ft 2 (l + -)-l = - 18 

ro 1 — c z 

As oj m > 0, so c 2 is restricted by the relation 

1 - b 2 < c 2 < 1 (19) 

However, in the plane (pd, p m ) the fixed points lie on a straight line with positive slope (-7 — !)• The 
phase portrait has been presented in Figure (I). In particular, if the DM is chosen in the form of Dust 
(i.e, Cold DM), then we have b 2 = 1 — c 2 , ro = ^ and hence the coincidence problem will have a 
solution by proper choice of the coupling parameters. 
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Figure I: The figure shows the phase portrait around the line of fixed points, where for the solid and 
dashed curves the initial values of (p m ,Pd) are (1, 0.1) and (0.8, 0.2) respectively. For both the curves 
we have taken b = 0.85, u m = 0.4 and c = 0.6. 



Case II: Q = 3AHp d 

This choice of interaction term appears to be phenomenological but it is compatible with observa- 
tions [32] like Snla, CMB, lagre scale structure, H(z), age constraints and recently in galactic clusters 
observations. Here A is a small dimensionless positive parameter and factor '3' is chosen for simplicity 
of calculations. 

In this case, the energy conservation equations have the explicit form 
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p! m = 3p rn [ v - (1 +(jJ m )] 

As before they form an autonomous system for constant uj m in the (pd, p m ) plane and the fixed 
points are characterized by 

w m = ^ - 1, r = ji-l 

As co m > 0, so the parameters are restricted by the relation 

I^<c 2 <l 

where the equality sign holds for cold DM. Similar to the previous case we have a line of fixed 
points with positive slope — in (pd, p m ) plane as shown in Figure (II). 



■* / 




/ V 


/ 


'/// 




\ 


\' 


\ 


vj v. _J ^ J ' 




/ 


/ / 


/ 






\ 


\ 


\ 






/ 


/ / 


/ 


fi 




\ 


\ 






f / 


/ 


/ / 


/ 






\ 


\ 






/ / 


/ 


/ / 


/ 




\ 


N 






/ / 


/ 


/ / 


/ 


ft \ 


\ 


\ 








' / 


/ 


/ / 


/ 


ft \ 


\ 


\ 








> / 


/ 


/ / 


/ 


i i 

t \ 


\ 


V 








* / 


/ 


/ / 




\ 


V 








■> / 


/ 


/ / 


\ 


t \ 


\ 


V 


- 






* / 


/ 


/ / 




\ \ 


N 


* 








' / 


/ 


/ / 


t 


\ \ 












s / 


/ 


/ // 


t 


\ \ 












s / 


/ 


/ 




\ V 


- 










' / 


/ 


/ ft 


\ 


\ V 












s / 


/ 


// 


\ 


\ ^ 












' / 


/ 


I \ 


\ 


s 












/ / 


/ 


t \ 


\ 














s / 


/ 


/ ^ 


\ 














' / 


/ 


// \ 


N 














•> / 


















' / 




K 


- 














' / 


\ 


\ V 
















s / 




\ V 
















* / 




\ - 
















* / 


















✓ / \ 


















n 
























































r - 1 - "«~ " " " " . . 


























— 1 




— i 










-_i -_i -_r _ 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

K 

i 

Figure II: The figure gives an idea of the phase portrait around the line of critical points where for 
the solid and dashed curves the initial values of (pmi Pd) are (0.8, 0.2) and (1, 0.1) respectively. For 
both the curves we have taken A = 0.6, uj m = 0.5 and c = 0.6. 

In the above discussion for both the cases the self-autonomous system can be written as pi = 4>{p) 
with p = (p m , pd)- To determine the stability criteria around the fixed points we expand the p around 
po (critical point), i.e, p = po + x, where the perturbation of the variables x satisfies (upto first order) 
xf = Mx. The linearized perturbation matrix M characterizes the fixed points and the stability criteria 
is determined by TrM < and detM > 0. However, in the present problem we have detM = (in 
both the cases), so the fixed points are degenerate equilibrium points. Further, the eigenvalues of the 
linearized matrix M are and Tr M, so the fixed points are non-hyperbolic in nature and it is not 
possible to obtain any conclusive information about the stability from linearization- one needs Normal 
forms calculations [33] or numerical experimentations. 
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4 A general prescription for deceleration parameter: 

Introducing the density parameters Q m = jff? and Q d = and the deceleration parameter q = 
— (1 + -§t), the field equations (2) and (3) can be written as 

n m + n d = i + n K (20) 

and 

1 3 

q = -(1 + Q K ) + -(fl d u d + tt m u m ) (21) 

with $7 K = 

Now solving for Q m and fi^ we obtain 

n m = [2q - (1 + 3w d )(l + fi«)]/3(w m - w d ) (22) 

and 

O d = [(1 + w ro )(l + n K ) - 2q]/3(uj m - co d ) (23) 
Hence the positivity of the density parameters restrict the deceleration parameter as 

1 -^(i + n K )< q < 1 -±p^(i + n K ) (24) 

which for the ACDM model (i.e, oo m = 0, oo d = — 1, k = ) simplifies to — 1 < q < \ and describes the 
whole evolution from dust era to ACDM. On the otherhand, during DE dominance (i.e, \ < U d < 1, 
< Vt m < |, with U K = 0) we have from (21) 

1 + 3u> d 1 3 , . 

-^-^<q< - + -(u; m + u; d ) (25) 

which on the phantom barrier (i.e, at the extreme DE era) simplifies to— 1 < q < — \. 

If we consider the system of two dark fluids as a single fluid with energy density pt = p m + Pd, then 

the effective equation of state of the single fluid is 

"* = i + n K (26) 

which does not depend on the interaction term. Hence the deceleration parameter simplifies to 

q = ±(l + 3u H )(l + n K ) (27) 

which is the usual definition of the deceleration parameter for a single fluid. 

Further, using the effective equation of state parameters from equations (6), the continuity equa- 
tions (4) and (5) can be integrated to give 

Pm = ^a(u) and p d = (28) 
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where u is a dimensionless variable defined as [34] 



a 



u=—, (29) 



a.. 



a* = a(i*) is the value of the scale factor at the instant t* when energy densities of both the dark 

components are identical (i.e, p m (i*) = p d (t*) = po(s&y)). 

Here, 

a(u) = exp[-3 H^^-dx] and (3{u) = exp[-3 f^Mrfx] (30) 

Ji x Ji x 



For flat space, the density parameters can be written as 

^ a(u) , ^ /3(u) . , 

On = n ^ n and n d = ^ 31 

and the deceleration parameter takes the form 

a(u) _ 1 | 3 Mu)Vm+P(u)u> d 



Now differentiating (32) with respect to u we get (after some simplification) 

P~ = ^(^m^J^ + ^d^) + 7rtt m ttd[-(Um ~ U d ) 2 + 7^7 (— + -)K - W d )] (33) 

2 au era 2u 3ii p d p m 



As w m , the equation of state parameter for DM is always positive and decreases (or remains con- 
stant) with the evolution of the universe, so < O.On the otherhand, the DE equation of state 
parameter oj d is always negative and can have any value (i.e, +ve or -ve). Thus 4^ may change 
sign more than once during the evolution of the universe and hence it is possible to have more than 
one real solution of q(u)=0. Therefore, it is possible to have more than one transition from decelera- 
tion to acceleration and vice versa in course of the evolution of the universe. However, in absence of 
interaction (i.e, Q = 0) if < then q is a strictly decreasing function of u [34] and q(u) = has 
exactly one real solution. Hence the interaction term has a significant effect in the transition from 
decelerating phase to accelerating phase or otherwise. 



5 Cosmological models: Some examples 

For flat model, from the equation (21), q = gives 

n m u m + M d uj d + i = i.e, lot + ^ = (34) 
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Also from the above equation we have 



Ud = 1 _o ( 35 ) 

we shall now consider the following cosmological models and examine whether more than one transi- 
tion is possible or not. 



5.1 Modified chaplygin gas: 

The equation of state for modified chaplygin gas (MCG) is given by [35] 



V = 19-^z (36) 



where < 7 < 1 and B, n > 0. Now, if MCG is considered as the combined DM and DE then the 
equation of state parameter is given by 

wr = - = 7 - -— 1 (37) 



so from equation (34) we have 



P n+l = -^T (38) 



which has only one feasible solution. Thus, there is only one transition from deceleration to acceler- 
ation as predicted by observations in recent past. Further, from the combined energy conservation 
equation, i.e, from p + 3H(p + p) = 0, for MCG model, we have on integration 

P n+1 = ^llC + £] 

where, p =3(7+1) (n+1) and C is the constant of integration. Hence the scale factor at the 
transition point is given by 



a=l gfr + j) 1* 

provided the constant of integration 'C is restricted by the relation C < . 



(39) 



Moreover, if MCG is chosen as the DE component then using the equation (35) we have 

B 

7 + — 



P n+1 = (40) 



Thus in any case we have only one transition from deceleration to acceleration to match with the 
present observation. 
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5.2 Holographic DE model: 



If the DE is chosen in the form of holographic DE with equation of state parameter as equation (13) 
then from equation (35) we have 

i n , VH^\ , Q _ ^m(i-fld)+| 
3 V 1 "I" c J ~T~ 3Hp d ~ Q d 

This simplifies to the cubic equation (choosing x = x/H^ ) 

2 1 1 

— x 3 + (-+a; m )x 2 + (6 2 -a; m --) = (41) 



when Q=3b 2 H(p m + 



while we have the cubic equation 



Y c xZ + ( A + \ + w «) x2 - K + ^) = ° ( 42 ) 



when Q = 3XHpd. 

It should be noted here that both the cubic equations will have exactly one positive root (provided 
b 2 < oj m + |) irrespective of whether oj m = or not. So, in this case also we see that only one 
transition is possible. 



ABp 



2a- 1 



5.3 The phenomenological model of the form p = — p — A a -i + p - 

This is a phenomenological model in which singularities can appear in the past or in future, depending 
on the choice of the parameters 'A' and 'B'. The following table shows the restrictions on the param- 
eters of the model and the physical quantities that diverge, while the other quantities remain finite 
[36, 37]. 

Table I 



Singularity 


Divergences of the physical parameters 


Restrictions on the model parameters 


Type I 


a — > oo, p — > oo, p — > oo 


| < a < 1, V A, B 


Type II 


p — > oo 


a < or -§ > 1 


Type III 


p — > oo, p — > oo 


a > 1 V A, B 


Type IV 


p — > oo 


< a < 5, V A, B 



Note that for a = 1, we have linear equation of state p = up with ui = — 1 — ^P B , which is 
basically a constant and it is not of much interest. However, if B < then strong energy condition, 
weak energy condition and null energy condition are satisfied but the dominant energy condition is 
violated here. Further, if a > 1 then we see that 

2a -1 



p — > —p — Ap"~ - as p 
and p — > —p — Bp a as p 







oo 



Thus u 



1 — 0(or — 1 + 0) as p — > and A > 0(or < 0), while uj — > +oo(— oo) when p — > oo 



and B < 0(or > 0). 



If the DE component is chosen as the above fluid then the equation of state parameter takes the 
form 

Ud =-i- ABp r (43) 

Ap^+B K ' 
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and from the equation (35) we have the quadratic equation 

ABx 2 - A5x - B5 = 



(44) 



where x 



-a—1 



and 5 



Q. m (l+UJ m )— I 



This quadratic equation has two positive real roots provided 
I. B < 



II. O m < minimum 
and 



III. 



4B 2 



< minimum[i, (1 + w TO )] 



4S 2 



Thus under the above restrictions it is possible to have two transitions for q(u) as shown in Figure 
(III) Further, as in the present model uj m > 0, restriction (III) implies B 2 < ^. Moreover, if the DM 

and B 2 is chosen to be A/8, then the restriction (II) becomes 



is chosen in the form of dust, i.e, u) Ti 

|, which is well within the observed estimation of DM component. 



Q. m < minimum^, w 



5r 




Figure III: These graphs show the variation of the deceleration parameter (q) with the dark energy 
density (p^) for the phenomenological model described in 5.3. 



6 Discussion and concluding remarks: 

The present work is a model of interacting dark matter and dark energy where both the matter com- 
ponents are in the form of perfect fluid having barotropic equation of state. If due to interaction 
there is a energy flow from DE to DM then the effective equation of state parameters for both the 
fluids change significantly- DM has a tendency to become exotic while the effective DE may be still 
in quintessence era even when we have phantom DE. The evolution equations for energy densities of 
the dark species from a self- autonomous system provided the equation of state parameter for DM is 
taken as constant. We have a degenerate line of critical points which are non- hyperbolic in nature 
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and hence the stability criteria can not be determined by linearization technique. 
In the subsequent section we have evaluated a general expression for the decelertion parameter in 
terms of the density parameters and the equation of state parameters-it does not depend on the 
interaction term explicitly. However, the evolution of the deceleration parameter depends on the in- 
teraction term as well as on the evolution of the state parameters. As from physical consideration, 
the DM equation of state parameter uj m should be positive but decreases with the evolution of while 
DE equation of state parameter uj^ is chosen as negative, so depending on the nature of the interac- 
tion term there may be more than one sign change in the expression 4^ in course of the evolution of 
the universe. Thus, it is natural to speculate that the universe will again be in decelerating phase 
after the present accelerating era. Finally, this speculation has been examined for some well known 
DE models and it is found that future decelerating era is possible only for the phenomenological model. 

Acknowledgements: Author SC is thankful to the UGC-DRS programme and author SP 
is thankful to The Council of Scientific and Industrial Research (CSIR), Govt, of India for Junior 
Research Fellowship. 
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